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Abstract 

A class of (l+l)-dimensional nonlinear boundary value problems (BVPs), modeling 
the process of melting and evaporation of solid materials, is studied by means of the 
classical Lie symmetry method. New definition of invariance in Lie's sense for BVP is 
(NJ \ presented and applied for the class of BVPs in question. 

O; 

1 Introduction 

> 

o; 

Boundary value problems of the Stefan type are widely used in mathematical modeling a 
huge number of processes, which arise in physics, biology, and industry pQ, [5], [11], [15], [T6] . 
Nevertheless these processes are very different from formal point of view, they have the common 
peculiarity, unknown moving boundaries. Movement of unknown boundaries is described by 
the Stefan boundary conditions. Boundary value problems with the Stefan conditions are the 
main object of this paper. It is well-known that exact solutions of BVPs of the Stefan type 
can be derived only in exceptional cases and the relevant list is very short at the present time 
(see P, [1], [6], [9], [T7] and papers cited therein). 

The main idea is to apply the classical Lie symmetry method ( [3], p2] [E]) to the case of 
the Stefan type BVPs. The Lie symmetry method is the very powerful tool for investigation 
of nonlinear partial differential equations (PDEs), notably, for constructing exact solutions. 
On the other hand, this method is efficient for solving standard BVPs (no moving bound- 
aries) in exceptional cases because the boundary conditions are usually not invariant under any 
\Q • transformations. 

Nevertheless the Stefan type problems are more complicated objects than BVPs with fixed 
boundaries, we have recently noted that the Lie symmetry method should be more applicable 
just for solving problems with moving boundaries [SJ. In fact, the structure of such boundaries 
may depend on invariant variable(s) and this allows us to reduce the given BVP to one of lower 
dimensionality. It should be noted that the rigorous definition of Lie invariance for BVPs and 
relevant examples are presented in [2UB]- However, this definition does not suit to many realistic 
BVPs of the Stefan type and should be modified and generalized. 

Let us consider the following class of the Stefan type BVPs used to describe melting and 
evaporation of materials in case that their surface is exposed to a powerful flux of energy 

l(x 2 (T 2) f)^ M (T 2 )f, (2) 

BT 

S l (t,x) = 0: X 1 (T 1 )-^ = p 1 (T 1 )L v V 1 -Q(t,T 1 ), V 1 = H(t,T 1 ), (3) 

BT BT 

S 2 (t,x) = 0: X 2 (T m )-^- = X l (T m )-^ + p 2 (T m )L m V 2 , T 1 = T 2 =T m} (4) 

x = +oo : T 2 = T , (5) 
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where T m , To are the known temperatures of melting and solid phase of material, respectively; 
Afc(Tfc), Ck(Tfc), and pk(Tk) are positive thermal conductivities, specific heat, and densities, re- 
spectively; L v , L m are latent heats of phase change per unit mass of liquid and solid phases, 
respectively; Q(t, T{) is a function that represents the energy flux being absorbed by the metal; 
H(t,Ti) is a given function; Sk(t,x) are the phase division boundary surfaces to be found; 
Vk(t,x) = ~^§f/^- are the phase division boundary velocities; Tk(t,x) are unknown temper- 
ature fields; index k — 1, 2 corresponds to the liquid and solid phases, respectively. 

Here Eqs. (JIJ) — (J2J) are basic and describe the heat transfer process in liquid and solid 
phases, the boundary conditions ([3]) present the evaporation dynamics on the surface Si, and 
the boundary conditions @ are the well-known Stefan conditions on the surface S2 dividing 
the liquid and solid phases. Assuming that the liquid phase thickness is considerably less than 
the solid phase thickness, one may use the Dirichlet condition ([5]). It should be stressed that 
we neglect the initial temperature distribution in solid phase and consider the process on the 
stage when two phases have already taken place. 

One may claim that formulae ([Q-flS]) present the class of BVPs with moving boundaries and 
take into account a number of different situations, which occur in the melting and evaporation 
processes. For example, setting Q(t,Ti) = const, V\ = const $(£) and H(t,Ti) = <3>(i)Ti, where 
$(t) is a correctly-specified function, one obtains the problem, which is the most typical [6], [9]. 
In the case of the process when surfaces are exposed to very powerful periodic laser impulses 
these functions take a complicated form jTU] . 

The paper is organized as follows. In section 2, we discuss the existing definition of Lie 
invariance for BVPs and present its generalization to a wider class of BVPs. In the next 
section, we apply the definition derived to BVP ([I])-©. In section 3, all possible Lie operators 
of the nonlinear system (GO)-©, which allow us to reduce the problem in question to one for 
ordinary differential equations, are found. Finally, we present conclusions in the last section. 



2 Definition of Lie symmetry for BVPs 

We start from the well-known definition of invariance of a BVP under the given infinitesimal 
operator presented in [3] . We restrict ourselves to the case when the basic equation of BVP is 
a two-dimensional evolution PDE of fcth-order (k > 2). In this case the relevant BVP may be 
formulated as follows: 

u t = F (x, u,u x ,..., ui k) ) , (t,i)eflc R 2 (6) 

s a (t, x) = : B a (t, x, u,u x ,..., uf~ l) ) = 0, a = 1, 2, . . . , p, (7) 

where F and B a are smooth functions in the corresponding domains, Q is a domain with 
smooth boundaries and s a (t,x) are smooth curves. Hereafter the subscripts t and x denote 
differentiation with respect to these variables, ui^ = §^j,j = 1,2, ... ,k. We assume that BVP 
©-(El) has a classical solution (in a usual sense). 
Consider the infinitesimal generator 

d d d 

x = e(t, x)-+ e(t, x)^+ r,(t, x, u)—, (8) 

(hereafter I; ,!; 1 and rj are known smooth functions), which defines a Lie symmetry acting on 
both (t, x, u)-space as well as on its projection to (t, x)-space. Let be the /cth-prolongation 
of the generator X calculated by the well-known prolongation formulae (see, e.g. [12]- |14j). 
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Definition 1 [3] The Lie symmetry X ftEj) is admitted by the boundary value problem (E|)-([7|) 
if and only if: 

(a) X^ (^F (x, u, u x , . . . , u x k ^ — Utj = when u satisfies (0|); 

(b) X(s a (t, x)) = when s a (t, x) — 0, a — 1,2, ... ,p; 

( c) X( k -V (B a (t, x, u, u x , . . . , u x k ~ l) ^j ) = when B a \ Sa{tjX)=0 = 0, a = 1, 2, . . . , p. 

One easily notes that Definition 1 can not be applied to BVP (CE])-©, because, firstly, there 
are two basic equations ([I])-© instead of one, secondly, the boundary condition ([5]) is defined 
on the non-regular manifold x = oo, thirdly, there are several boundary conditions on the 
moving surfaces Si(t, x) and S2(t,x). Thus, this definition should be extended to the wider 
class of BVPs. 

Consider a BVP for the m-component system of evolution equations with two independent 
(t, x) and m dependent u = (ui,U2, ■ ■ ■ , u m ) variables. Let us assume that the fcth-order (k > 2) 
basic equations of evolution type 

u\ = F i (x, u,u x ,..., u x k) ) , i = 1, 2, . . . , m, (9) 

are defined on a domain Q C R 2 and there are three types of boundary conditions, which can 
arise in applications: 

s a (t,x) = : B'i (t,x,u,u x ,...,u x k ^) =0, a = 1, . . . ,p, j = 1, . . . , n a , (10) 

((99 (99 \ 
t,x,u,...,u x k »\S b ,-^,-^\ =0, b = l,...,q,l = l,...,n b , (11) 

and 

7 c (t, x) = oo : T c (t, x, u,u x ,..., u ( x c) ) =0, c = 1, 2, . . . , r. (12) 

Here kj a < k and < k are the given numbers, s a (t, x) and 7 c (t, x) are the known functions, 
while the functions Sb{t,x) defining free boundary surfaces must be found. We assume that all 
functions arising in (I9~j)- ([T2"|) are sufficiently smooth so that a classical solution exists for this 
BVP. 

Consider the infinitesimal generator 

X = e(t, x)^ + £% x )± + r ,% x ,u)A. + ... + v -(t, x, u)-^, (13) 

which defines a Lie symmetry acting on both (t, x, u, >9)-space (here the notation S = 
(Si, 5*2, ... , S q ) is used) as well as on its projection to (t, x)-space and generates a Lie group of 
the point transformations 

t' = T(t,x,e), x' — X(t,x,e), u'j = Uj(t,x,u,e), S' b = Sb(t,x). (14) 

Definition 2 BVP (fPj)-(fl^) admits the one-parameter Lie group of transformations p^j j gen- 
erated by the infinitesimal operator H~3}) if and only if: 

(a) X^ (^F l (x, u, u x , . . . , u x k ^ — u\j = when the functions Ui, i = 1,2, ... ,m satisfy (G|) ; 

(b) X(s a (t, x)) = when s a (t, x) = 0, a = 1,2, . . . ,p; 
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(c) [B{ (t,x,u,u x ,...,u£ ja) )) = when Bi\ Sa(t>x)=0 = 0, a = 1,2, ...,p, j = 

(d) X^ ^B l b (t,x,u,...,u x klb) ,S b ,^,^ = when B l b \ Sb{t>x)=0 = 0, b = l,...,q, I = 
1, . . . ,n b ; 

(e) x *(lc( T ,y)) = whenY c {r,y) = 0, c=l,2,...,r; 

(f) Xi hc) (r* (t, y, u,u y ,..., 4 fcc) )) = when r*| 7 * (T)2/)=0 = 0, c = 1, 2, . . . , r, 

where X*, Yl and 7*(t, y) are operator ( E5p , the functions Y b and - ^ tx y respectively, expressed 

( t, if £>7c(f ' a:) ^ 0, 

ma t/ie new independent variables r = < ' x ^foU) g' an< ^ 1/ = 

^ 7 C (i, a;)' ^ 9x 

if T = ' 
i z/ ajcM ^ Q 

It should be noted that this Definition coincides with Definition 1 if m = 1, n a = 1 and 
there are no boundary conditions of the form (TTTi) and (|T2|) . On the other hand, the following 
example shows that the generalization of Definition 1 presented above is non-trivial. 

Example. Let us consider BVP (!9l)- (jT0{) . which includes also the boundary conditions 

x = oo : T c (u) = 0, c = 1, 2, . . . , r. 

Assume that system d§J) and conditions (TTOi) are invariant under the group of translations on 
the plane (t, x): 

t' = t + \xE, x' = x + A 2 £, it' = u, AiA 2 7^ 0, (15) 
so that the corresponding infinitesimal generator 

X = \ x d t + X 2 d x (16) 

satisfies items (a)-(c) from Definition 2. Moreover, according to Definition 2 the operator X of 
the form ffT6l) and the functions r c (ii) and j c (t,x) = x take the form 

X* = Ai<9 T - \ 2 y 2 d y 

r*(tt) = r c (tt) and 7*(t, y) = y, respectively. Now one easily checks that items (e) and (f) from 
Definition 2 are satisfied. Thus, BVP in question admits the Lie group of point transformations 
()15p generated by the infinitesimal generator ffTBP . 

However, the problem with non-regular manifold 7 c (t, x) = oo occurs if one generalizes 
Definition 1 in the standard way and formulates items (e) and (f) from Definition 2 like those 
(b) and (c). For example, this manifold can not be replaced by the regular one 7 c (t, x) — L = 0, 
where L — > oo because this leads to the requirement 

gi%(M)-% H = o. 

On the other hand, this requirement is not satisfied in the example presented above because 

lim X(x - L)\ L = lim A 2 = A 2 + 0. 

L— >oo L— >oo 

Thus, item (e) in Definition 2 can not be formulated like item (b). 

Hereafter Definition 2 will be applied to derive the invariance operators of BVP ([I])-©. 



4 



3 Invariance of BVP (JTj)— (Ej) under the Lie generators 

It should be noted that BVP (pQ)-([55 can be simplified if one applies the Goodman substi- 
tution 

u = MTi) = \ ci(C)pi(C) d(, v = 2 (T 2 ) = [ c 2 (£)p 2 (£) d£. (17) 



Substituting (I17p into JI])-© and making the relevant calculations, we arrive at the equivalent 
BVP of the form 

« = te(* (u >a;)' (18) 

(19) 



<9t 9a; \ <9x 

S'i(t,x) = 0: = pi(u)L v V 1 -q(t,u), V 1 = h(t,u), (20) 

S 2 (t,x) = : d 2 {y m ) — = dx(u m )— + p 2 (v m )L m V 2 , u = u m ,v = v m , (21) 
x = +oo : f = Voo, (22) 

To 

where u m = J Ci(C)pi(C) d(, v m = J c 2 (£)p 2 (£) d£, Woo = / e 2 (£) p 2 (£) d£; = 

^ifc^ " <W = ^(i))" = PiC^M), h{v) = P^ 2 \v)), q(t,u) = 

Q(t, 0^ 1 (m)), and h(t,u) = 0J~ («)) (here 0" 1 are inverse functions to fa, the functions 
di(u) and <i 2 (v) are strictly positive and v m ^ Voq). 

Now one sees that BVP (I18p - (l22p is based on the standard nonlinear heat equations (NHE). 
Lie symmetries of the non-coupled system (j!8p - (ll9p can be easily derived from paper [7], where 
reaction-diffusion systems of more general form have been investigated. The result is presented 
in Table 1 (note, we do not consider the case when this system is linear). Using Definition 2 
and Lie symmetries from Table 1 the following theorem can be established. 



Theorem 1 A nonlinear BVP of the form fufy - tzty) admits a Lie symmetry operator of system 
(flff p -jTIP p if and only if the operator in question up to the local transformations x — >■ x + x , t — > 
t + t Q (xQ e R, t E R) is equivalent either to the operator 

X 1 = d t + /id x , /jgR (23) 

or to 

X 2 = 2td t + xd x . (24) 
Moreover, the functions f(t,u), q(t,u) must have the correctly specified forms 

h(t,u) — h(u), q(t,u) — q(u) (25) 

and 

7 , \ h(u) „. . q(u) , , 

h(t,u) = -^J-, q(t,u) = ^l, (26) 

respectively. Here q(u) and h(u) are arbitrary smooth functions. 
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Proof. Firstly, we derive the group classification of the NHE system ()l8]) - (jl9p using the 
determining equations presented in [7]. If di(u) and d 2 (v) are arbitrary functions then those 
equations immediately give the three-dimensional maximal algebra of invariance (MAI) A = 
(d t , d x , 2td t + xd x ) called the principal algebra. There are five special cases when an extension 
of the principal algebra occurs and they are listed in Table 1. 

It should be stressed that each NHE system admitting four- or five- dimensional Lie algebra 
is reduced to one of those with diffusivities from Table 1 by the equivalence transformations 



t — > e t + t , x — > e\x + x , u — > e 2 u + u , v — > e 3 v + v , 



(27) 



where ^ (i — 0, . . . , 3), to, Xo, Uo, and vq are arbitrary parameters. It turns out that the 
class of BVPs ( fl~8]) -( l22"]) is also invariant under transformations f )27|) and this can be easily 
checked by direct calculations. 

Thus, on the next stage, we need to apply Definition 2 only to BVP ( TT8"j) -( l2"2"j) in six cases 
listed in Table 1. 

Table 1. Lie algebras of the NHE system flT8l)-flT9lL 



no 


di(u) 


d 2 (v) 


Basic operators of MAI 


1. 


V 


V 


A=(d t ,d x ,2td t + xd x ) 


2. 


V 


1 


A, vd v , f3(t, x)d v 


3. 


e u 


e° 


A, xd x + 23 u + 2d v 


4. 




v m 


A, xd x + 2d u + ±vd v 


5. 


u n 


v m 


A, xd x + \ud u + ^vd v 


6. 


4 

U 3 


4 

V 3 


A, xd x — |ui9 u — |t> d v , x 2 d x — 3xud u — 3xvd v 



Let us consider the first case when the functions d\{u) and d 2 (v) are arbitrary. In this case 
the most general form of the Lie symmetry generator is 



X = (Xi + 2X 3 t)d t + (A 2 + X 3 x)d x 



(28) 



(hereinafter A with indices are arbitrary constants). Obviously, operator (12"8"|) can be reduced 
either to the form (123]) (if A3 = and Ai 7^ 0) or to (|24|) (if A3 7^ 0) by the local transformations 

x — > x + x , t — > t + 1 (x g R, t e R) , 

which belong to (|2T|) . 

Note that the case A3 = Ai = leads to the invariance operator X = d x . However, this 
operator generates the Lie ansatz, which does not alow existing the moving boundary surfaces 
Sk(t, x), k = 1, 2, hence, we do not consider this non-physical case. 

Let us apply Definition 2 to prove the invariance of BVP (fT8"|) - (l2"2"|) under the Lie symmetry 
operator (12"3"|) . Applying the first prolongation of generator ^ n v " 
conditions (1211 . we obtain 



Xi [d 1 {u)—-p 1 {u)L v V x + q{t,u, 



X\ = Xi, to the boundary 
dq(t,u) 



M 



dt 



Xi [Vi-h(t,u] 



dh(t, u) 
dt~' 
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where the manifold M, — {iSi(t, x) — : g?i(u)|^ = p\(u)L v Vi — g(t, it)} and the manifold A/" = 
x) — : Vi — h(t,u)^. According to item (d) from Definition 2 we arrive at the restric- 
tion on the function q(t,u) and h(t,u) 

dq(t, u) , . 

J^_J_^ = ^g = g( M ), 

at 

The invariance of conditions fl2T|) are trivially fulfilled because they do not involve the variables 
t and x in explicit form. Finally, condition f )22|) satisfies items (e) and (f) from Definition 2 
because this immediately follows from Example. 

The invariance of BVP (fTBj) - ([2"2"j) under generator ([23]) is checked in a similar way. Note 
that the first prolongation of generator (12"4"|) has the more complicated form than one of (|23|) . 
namely 

XP = 2td t + xd x - 2u t d ut - 2v t d vt - u x d Ux - v x d Vx - 2S\d sl - 2S 2 d s? - S%i - S 2 x d sl (29) 

Since the boundary conditions (1201) involve the first-order derivatives (we remind that Vk = 
~^§t/^§x) we have calculated how operator (!29j) acts on these boundary conditions according 
to Definition 2 and obtained two ODEs of the form 

dq(t,u) .. , dh(t,u) ~. 

2tJi ^ + M ) = °> 2t m + M ) = °> 

which lead to the restrictions (|26p . 

Thus, the case of arbitrary functions d\{u) and d 2 (v) is completely examined. 

It turns out that any other Lie generator arising in Table 1 does not satisfy Definition 2 for 
BVP (fI5l) - (j2"2l . Let us consider, for instance, case 4 of Table 1. Here, the most general form 
of the Lie symmetry generator is 

X 3 = (Ai + 2X 3 t)d t + (A a + (A 3 + A^a;)^ + 2A 4 <9 U + -X A vd v , (30) 

m 

where we assume A4 7^ 0, otherwise X 3 takes form f l28|) . 

In order to be invariant under generator fl30|) . the second condition from f l2T|) must satisfy 
item (d) from Definition 2 

X 3 (u - u m )| p = 0, 

where V = {^(t, a;) = : w = w m }, however, this leads to the restriction A4 = 0. Thus, the 
contradiction is obtained and we conclude that generator fl30|) is not Lie symmetry of any BVP 
f n8|) -f l22|) with diffusivities from case 4 of Table 1. 

In a similar way, cases 2,3, 5 and 6 listed in Table 1 have been examined and contradictions 
with boundary value conditions established. 

The proof is now complete. ■ 

Remark 1 The group classification of the NHE system / T7£]) -/ f7P)) presented in Table 1 takes 
into account the trivial discrete transformations u — > v , v — > u admitted by this system. Never- 
theless, the class of BVPs (T72j)-([II)) is not invariant under these transformations (see conditions 
fiEB) ) it does not affect Theorem 1. 
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Using Theorem 1 one can reduce the class of BVPs fTT8]) - f[22]) to two classes of BVPs for 
second-order ODEs. In fact, substituting ansatze 

u = «(£), v = v(£), Sk = Sk(£), £ = x-ixb, (31) 

X 

u = u(u), v = v(uj), S k = S k (u)), uj = —=, (32) 

yt 

generated by operators fT23|) and f )24|) . respectively, and taking into account restrictions f |25|) 
and ( 12"B"|) the following consequences are obtained. 

Corollary 1 Ansatz l[3l\) reduces a nonlinear BVP of the form with the coefficient 



restrictions l[25*\) to the following BVP for the second- order ODEs 

d ( , , .du\ du „ . . 

^(* (!,) 5?) + ^ = °- 0<f<i ' < 33) 

d ( , / n dv\ dv . 

^(^h' 1 ^ ' (>e - (34) 

£ = 0: d l {u)—=p 1 (u)L v ii-q{u),iJ, = h{u) 1 (35) 

£ = 5 : d 2 (v m )^ = di(u m )^ + p 2 (v m )L m p, u = u m , v = v m , (36) 
at, at, 

£ = +oo : v = foo, (37) 
where 5 and fi are parameters to be found. 



Corollary 2 Ansatz reduces a nonlinear BVP of the form / f7^) -/ f^j) wtt i/ie coefficient 



restrictions l[26}) to the following BVP for the second-order ODEs 

d { . , s du\ uj du n , 00 . 

— d^u)— + -— = 0, U! < uj < uj 2 , (38) 
duo \ duo J 2 du 

d f , / \ dv \ uj dv , , 

^(^) + 2^ = °- W>W2 ' < 39 > 

. . G?1t i i CiJi . . UJ\ , . . . 

w = : "iW^j = Pi\ u )L v — - q[u), — = h{u), (40) 

, , \dv , , \ du „ uj 2 

uj = u 2 : d2{v m )— = di{u m )-^ + P2\y m )L m — , u = u m , v = v m , (41) 

cj = +oo : v = Voo, (42) 
where to\ and uj 2 are parameters to be found. 



It should be noted that each BVP fl33l - fl37j) with the correctly specified functions q(u) and 
h{u) arising in fl35|) can exactly be solved and an exact solution can be presented at least in an 
implicit form (see for details j6]). Some BVPs fl38|) - fj42|) with the correctly specified coefficients 
arising in the problem were exactly solved in particular cases [8]. 

4 Conclusions 

In this paper, the class of (1+1) -dimensional nonlinear boundary value problems (pQ)-([S]), 
modeling the process of melting and evaporation of metals, is studied by means of the classical 
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Lie symmetry method. New definition of invariance in Lie's sense is presented. This definition 
is applicable to wide classes of BVPs, including those with several basic equations, with moving 
boundaries, and with boundary conditions on non-regular manifolds. 

Theorem 1 giving all possible operators of Lie's invariance for the class of BVPs ([I])-© is 
proved and the corresponding corollaries, which allow us to reduce the problems in question to 
those for ODE systems, are also obtained. 

Finally, we note that Theorem 2 from the recent paper [8] follows as a particular case from 
Theorem 1, nevertheless Definition 2 was not used in [8] because the direct application of Lie 
point transformations was applied therein. This means that definition of Lie's invariance for 
BVPs can be derived in different ways. We are going to discuss this in details in a forthcoming 
paper. 
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